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Abstract A geometric approach is used to study the Abel first order differential 
^ ' equation of the first kind. The approach is based on the recently developed theory 

, of quasi-Lie systems which allows us to characterise some particular examples 

of integrable Abel equations. Second order Abel equations will be discussed and 
the inverse problem of the Lagrangian dynamics is analysed: the existence of two 
alternative Lagrangian formulations is proved, both Lagrangians being of a non- 
natural class. The study is carried out by means of the Darboux polynomials and 
. Jacobi multipliers. 
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1 Introduction: Abel differential equations 



> 

^SJ ■ Nonlinear differential equations play a relevant role in the study of evolution in 

OA ' terms of an evolution parameter, what has motivated physicists' interest during 

, the last forty years. The lack of a general procedure for determining their gen- 

eral solutions justifies the analysis of particular instances where solutions can be 
investigated via algebraic or geometric methods. An important example are Lie 
systems, which admit a superposition rule giving the general solution in terms of 
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a finite set of particular solutions. These systems describe the integral curves of t- 
dependent vector fields which are t-dependent linear combinations of vector fields 
closing on a finite-dimensional real Lie algebra (see e.g. [T] and [5] for a modern 
presentation of the theory of Lie systems _3 , 4_ ) . 

An interesting instance of Lie system is the Riccati equation, a generalisation of 
the inhomogeneous linear equation. One of its main characteristics is the existence 
of an action of the group of curves in SL{2, R) on the set of Riccati equations, 
what can be used for reduction of a given Riccati equation into a simpler one. This 
is for instance the case when a particular solution of a given equation is known 
[5] . Moreover, this property can be shown to be common to all Lie systems [6J . 

A generalisation of Lie systems has recently been proposed 7,8J allowing us 
to use some of the techniques applicable to Lie systems to this more general class 
of systems. This paper deals with a specific example, the Abel equation, with 
applications in many different branches of physics. 

There are several differential equations called Abel equation. The first-order 
Abel equation of the first kind 

± = Ao{t) + Aiit)x + A2{t)x^ + A3(^)a;^ (1) 

is a generalisation of the Riccati equation (the particular case A3 = 0) introduced 
by Abel within the theory of elliptic functions, while the first-order Abel equation 
of the second kind is of the form l9lllOI : 



{V + f[t)) y = Bo{t) + Bi{t)y + B2{t)y'' + B3{t)y\ (2) 
which reduces into one of the first kind by means of the transformation 

x = {y + f{t)r\ (3) 

with coefficients 

Ao = -B3, Ai=W3f -B2, 
A2 = -f -SBsf + 2fB2- Bi, 
A3 = fB3-f B2 + /B1-B0. 

Therefore, we can restrict ourselves to study the Abel equations of the first kind. 
We can also consider the generalised first-order Abel equation [llj 



J2^k{t)x'', n>3, neN, (4) 



X = 

k=0 

or the second-order Abel equation [12] : 

d^x . 2 f dx\ ^ 

The n-order Abel equation is a linear combination of the different members of a 
hierarchy of j-order Abel equations with functions pj = pj (t) as coefficients [12] : 

(poBa + PiDa"^ + ■ ■ ■ + Pu-iDa + Pn)x + p„+i = . (6) 

Here, Ha denotes the differential operator = d/dt + x^ , and the action of 
Dyi leads to a family of differential equations whose first members are given by a 
sequence D^a; = 0, j = 0, 1, 2, . . .. 
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Abel equations appear in the reduction of order of many second- and higher- 
order equations, and hence are frequently found in the modeling of real problems 
in many areas, for instance, Emden equation, y" + {2/x)y' = 0, Emden-Fowler 
equation, y" = A y"^ {y')\ the Van der Pol equation y" - e(l - y'^)y' + ay = 0, 
the Duffing equation, y" + ay + by^ =0, and many other equations. They also 
play a relevant role in the study of quadratic systems in the plane [13] and the 
center-focus problem [14] . 

The Lienard equation x" + f{x)x' + g{x) = was studied in jl5| . Defining 
^(x) = x' , it may be written -|- f{x)(, + g{x) = 0, which is an Abel equation 
of the second kind, related by u = l/y (see ([3])) to the Abel equation of the 
first kind: u' = f{u)u^ + g[u)u^. In particular, the one studied by Bougoffa [10] . 
5,5,' — £, = 'P{x) reduces to u' + + <P{x) = 0. For instance, the autonomous 
dissipative Milne-Pinney equation can be transformed into a second kind Abel 
equation: 

which is related to the first kind Abel equation x = x'^ — j^x^ ■ 

The Abel equation was used by Majorana in the study of Thomas-Fermi equa- 
tion: 

,,3/2 

y =—2^ (7) 

which is a particular case of Emden-Fowler equation, and it is also used in the 
study of cosmological models [16] . The general solution of the cosmological Einstein- 
Friedmann equations for the universe filled with a scalar field for a given potential 
can be expressed via the general solution of the Abel equation of the first kind 
[17] . Another more recent application can be found in [18j . 



2 Lie systems and superposition rules 



There are (systems of n) first-order differential equations (to be called Lie systems) 
that admit a superposition rule, i.e. a function 4) : M"("^+i) R" such that the 
general solution can be written as a; = (j)[x(^y^, . . . , a;^^), fc), the simplest example 
being the inhomogeneous linear equation 

x = co{t) + ci{t)x, (8) 

for which there is a superposition rule involving two generic particular solutions: 

x{t) = ct>{xi{t),X2{t); k) = xi{t) + k {x2{t) - xi{t)). (9) 

Moreover, it admits a solution in terms of two quadratures: 



c(t) = exp 



I Cl{t')dt' 



c 



exp 



Cl{t")dt" C2{t')dt' 



Another well known example is the Riccati equation, i.e. the following nonlinear 
differential equation: 



X = Co(t) -I- Cl{t) X + C2(t) X 



(10) 
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This type of equations does not admit integration by quadratures in the general 
case. Their solutions are the integral curves of the t-dependent vector field 

Xt = Co{t)Xo + Cl{t)Xl + C2{t)X2, 

which is a linear combination with t-dependent coefficients of the vector fields 

Xq = dx, Xi = xdx , X2 = dx , 

that close on a three-dimensional real Lie algebra, with defining relations 

[Xq, Xi] = Xq , [Xq, X2] = 2X1 , [Xi, X2] = X2 , 

isomorphic to the s[(2,R) Lie algebra. Therefore it is a Lie system admitting a 
superposition rule which turns out to be [19j 

(f,(xi,X2,X3;k) = — -j ' ' (11) 

fc (^3 - a;2) + [XI - X3) 

The vector fields Xq, Xi and X2 are a basis of the fundamental vector fields 
relative to the action of the group SL{2, R) on R = R U {00}: 



${A,x) 



( ax + P J r r / 1 

-, li x {— 0/7, 00}, / „s 

^^ + ^ withA=f"^)e5L(2,I 



if X = —6/"f, 



A Riccati equation can be seen as a curve in R'^ and each element A{t) of the group 
Q of smooth SL ( 2, R)- valued curves, Map(R, 5L(2,R)), transforms every curve 
x{t) in R into a new one x{t) = <l>{A{t),x(t)) satisfying a new Riccati equation 
with coefficients C2, ci, cq: 

C2 = C2 — 57 ci + 7^ Co + 7(5 — 57 , 

ci = — 2 C2 + {aS + l3"f) ci — 2 07 cq + 5a — a5 + f^-y — 7/? , 
Co = P'^ C2 — afi c\ + Q'^ Co + a/3 — fia . 

These expressions define an affine action of the group Q on the set of Riccati equa- 
tions. An appropriate choice for the curve, A{t), transforms the original Riccati 
equation into a simpler one, for instance an integrable by quadratures one. The 
conditions for the existence of such a transformation provide us with integrability 
conditions . Moreover, similar properties hold for any other Lie system. 



3 Geometric approach to the Abel equation of the first kind 

Each t-dependent vector field whose integral curves are the solutions of an Abel 
equation is a linear combination with t-dependent coefficients of the vector fields 
in the linear space VAbei(R) defined by 

VAbel(R) = (^dx,xdx,x'^dx,x^dx^ ■ 

Such a f-dependent vector field is not, in general, a Lie system. Indeed, consider a 
i-dependent vector field X(t,x) = x'^dx + h(t)x^dx, with a non-constant function 
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h. If X{t,x) determines a Lie system, there must exist a finite-dimensional Lie 
algebra Vo{R) including the vector fields x'^dx and x^da: such that X e Vo{C°°{R)). 
Nevertheless, this is impossible as such a Lie algebra should contain the successive 
Lie brackets of x^dx and x'^dx , which span an infinite family of linearly independent 
vector fields (note that [x'^ dx , x"' dx] = (n - 2)x"+'^dx)- 

As pointed out in a recent paper [7] , we can however look for a linear subspace 
W{R) C yAbei(]R) such that W(R) is a Lie algebra satisfying that [M^(R), FAbeiW] C 
VAbei(IK). In such a case the flows of time-dependent vector fields with values in 
14^ (R) leave invariant the linear space VAbei(R) and starting with a given Abel 
equation we obtain a new Abel equation from each flow. Maybe, by an appropri- 
ate flow, the transformed Abel equation becomes a Lie system. In this case the 
initial Abel equations is called a quasi-Lie system. 

Our aim is, given Y = f(x) dx G j£(K), to determine / such that [Y, VAbei(R)] C 
VAbet{R)- In particular, 

[Y,dx] = -f'{x)dx e ^AbelW- 

Hence, —f'{x) = C3X^ +C2X^ +cix+co and f{x) = —^x'^ — ^x^ — ^x'^+cqx+c-i. 
Using now the condition [Y,^^^;] = {fix) — xf'{x))dx G VAbeKl^), we obtain that 
C3 = 0, and \Y,x'^dx] = i2xf{x) -x^f'{x))dx G VAbei(K) yields that C2 = 0. 
Moreover, as Y must satisfy that \Y,x^ dx] = {3x'^f{x) - x^f'{x))dx £ VkbeKK), 
we also see that ci = 0. Consequently, every symmetry group of VAbei(K) admits 
a Lie algebra of fundamental vector fields Wa contained in the Lie algebra 

W^Abel = {dx,xdx) ■ 

We find in this way the so-called structure invariance group pP , which turns out 
to be the affine group in one dimension. 

Correspondingly, the set of first-order Abel equations of the first kind is invari- 
ant under all the following transformations [22] 

x{t) = a{t)x{t) + P{t), a{t)=^0, (12) 

Under such a transformation the given Abel equation ^ becomes a new Abel 
equation x = Ao{t) + Ai{t) x + A2{t) -|- A3{t) x^ , with 

A^it) = A3it)a^it), 

A2(t) =Q(t)(3^3(i)/3(^)+^2(t)), 

Ai{t) = 3A3{t)P^{t) + 2A2{t)l3{t) + Ai{t) - a{t)a-\t), ^^^^ 
Aoit) = a-\t) (A3{t)P\t) + yl2(i)/3'(i) + ^i(i)/3(t) + ^o(i) - /?(*)) • 

This is a very useful property: the action of the group of curves in the affine 
group produces orbits, i.e. equivalence classes of Abel equations, with the same 
properties of integrability or existence of superposition rules. For instance, Abel 
equations of an orbit containing an integrable by quadratures Riccati equation are 
also integrable by quadratures. The orbits are characterized by different values of 
invariant functions. 

One must determine the Lie algebras contained in VAbei reachable by the set of 
transformations (the coefficient of X3 cannot be transformed to be zero, see first 
equation in psp l. It is possible to check that the only subalgebra of dimension 
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three is the one of Riccati equation which is not reachable by the considered set 
of transformations. On the other hand, there is a one-parameter family of two- 
dimensional subalgebras of VAbei reachable from a proper Abel equation, those 
generated by 

{-2fi'^ +3lJ.X^ +x'^)da:, {fi + x)d^. (14) 

These are integrable by two quadratures. For instance the case fi = corresponds 
to Bernoulli equation. Finally, the algebras generated by vector fields 

{co + ClX + C2X^ + x^)dx, (15) 

lead to separable differential equations, therefore integrable by one quadrature. 



4 Canonical forms of Abel equations of the first kind 

Starting from an Abel equation of the first kind with ^ 0, we can define a 
transformation x = x + ^A2(t)/ A^lt), which reduces the original equation to a 

new one x = Ao{t) + Ai{t) x + A3 {t) x^ , where 

J , A1A2 2 Ai Id fA2\ J . 1^1 7 ,^ 7 

In fact, a translation x = x + /3 transforms the given equation into: 
S = /3 + A3{x - af + A2(x - of + Ai{x - a) + Aq, 

and then 

1 = A3X^ + {A2 - 3/3^3 )5^ + (3/J^ ^3 - 2/3^2 + Al )a; + /3 - A3 + A2 - + Ao . 

Choosing /3 = we find that A3 = A3, A2 = and Ai and Aq are given 

3A3(t) 

by the above given expressions. When Ao = 0, we obtain a first possibility for the 
canonical form 

X = Ai{t) X -\- Az{t) x^, 

which is a Bernoulli equation for n = 3, and therefore solvable by the change of 
variable u = 1/x^, which leads to the inhomogeneous linear equation 

u + 2Aiu + A3 = . 

On the contrary, if Ao_^ 0, under the transformation x = xAq, the equation 
becomes x' = A3{t)x^ + Ai(a;) y + 1, and we obtain as a canonical form 

X = 1 + Ai{t) X + A3{t) x^ . (16) 

An important instance is the case Aq = Ai = 0: the canonical forms are such that 

lAi - _ 2 Ai 1 d /A2\ 

SM' ^°-27A| + 3dt (.A^J- 

Liouville proved that if #3 and ^5 are given by 

#3 = A2 A's - A'2 A3 + 3 Ao Ai - Al A2 A3 + |Ai 
^5 = ^3 ^'3 - 3 (A?, + iAi - Al A3) 
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then the quotient ^3/^5 is an invariant under the structure invariance group. The 
first canonical form corresponds to ^3 = 0. Any two equations with ^3 = are in 
the same orbit, i.e. are related by the structure invariance group. 

For equations with the second canonical form there is a new invariant, and 
they are in the same orbit if and only if the value of the invariant is the same for 
both equations. 

There exist a two-parameter structure invariance group for this type of equa- 
tions in the case of first canonical form and a one-parameter one in the cases of the 
second form provided that an invariant has a constant value different from zero. 



5 Second order Abel equation 

Lie's recipe for order reduction of second-order linear differential equations leads 
to some Riccati equations. For instance y = under the change of variable y = 
becomes Da; = with D = d/dt + x , where x = z. Similarly y = becomes under 
the same change of variable D^x = 0. On the other hand, when comparing Abel 
equation with Riccati equation one sees that a possible generalization to second- 
order equation of Abel equation is 0\x = 0, where D^ = d/dt + x'^ [12]. The 
second-order Abel equation ([3]) can be presented as a system of two first-order 
equations 

dx 

dv 2 5 

— = —4a; V — X , 
dt 

corresponding to the following vector field on the velocity phase space 

r = vd:^+FAdv, Fa = -'ix^v-x^. 
Such a nonlinear Abel equation of second-order can be derived from the Lagrangian 

LAix,v) = ^ . (17) 

The corresponding conserved energy function is given by 

(3^ + a;3) 
(w + a;3)3- 

To be remarked that there exists an alternative Lagrangian given by: 

LA(x,v) = [?,v + x'^fl^. (18) 

Such Lagrangians can be obtained by means of the Jacobi last multiplier theory 
and Darboux polynomials for polynomial vector fields. More specifically, Darboux 
polynomials for a polynomial vector field X are polynomials I? such that XD = 
/P. The function / is said to be the cofactor corresponding to such a Darboux 
polynomial and the pair (/, I?) is called a Darboux pair. 

On the other side, given a vector field X in an oriented manifold (M, ^2), a 
function R such that Ri{X)f2 is closed is said to be a Jacobi multiplier (JM) for X 
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[23| . Recall that the divergence of the vector field X (with respect to the volume 
form Q) is defined by the relation 

Cx^ = (divX) r?. 

Then, i? is a Jacobi multiplier if and only if i?X is a divergenceless vector field 
and therefore using 

Crx^ = {divRX) n = \X{R) + i?divX] i? = , 

we see that i? is a last multiplier for X if and only if 

X(i?) + i?divX = 0. (19) 

because, for any function /, 

^(/ ^) + / div X = (X/) R + f{X{R) + R div X). 

The remarkable point is that if Di , . . . , 2?fe are Darboux polynomials with cor- 
responding cofactors fi, with i = 1, . . . ,k, one can look for multiplier factors of 
the form 

k 

R=Y[Vi\ (20) 

i=l 

and then 

X{R) _ ^ X(PO 



1=1 1=1 

Consequently, if the coefficients vt can be chosen such that 

k 



Y,V^f^ = -AiYX, (21) 



^t^i fi = -divX, 



we arrive to 

X{R) 
R 

i=l 

what implies that i? is a Jacobi last multiplier for X. 

Finally, one can prove (see e.g. ^23^) that if i? is a Jacobi multiplier for a 
vector field which corresponds to a second-order differential equation, there is 
an essentially unique Lagrangian L (up to addition of a gauge term) such that 
R = d'^L/dv^. 

Actually, Helmholtz analysed the set of conditions that a multiplier matrix 
gij(x,x) must satisfy in order for a given system of second-order equations 

= F'^ (x, x) , j 

or the corresponding system 

( dx' 
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which provides us the integral curves for the vector field 

X = v'd,.+F\x,v)d,., (22) 

when written of the form 

gi-j{x^ -F\x,x)) = 0, = 1,2, . . . ,n, 

to be the set of Euler— Lagrange equations for a certain Lagrangian L. 

Each matrix solution gij can be identified with the Hessian matrix of L, gij = 
dL/dv^dv-' , and then L can be obtained by direct integration of the gij functions. 

The two first conditions just impose regularity and symmetry of the matrix 
gij ; the two other equations introduce relations between the derivatives of gij and 
the derivatives of the functions F^ . The fourth set of conditions for gij is 

Xigzj) = gtkA^ j + gjkA^ i , = -^^d^jF' . 

When the system is one-dimensional we have i = j = k = \ and then the three first 
set of conditions become trivial and the fourth one reduces to one single P.D.E. 

, dF do ^ dq OF 
X{g) + g—=v^+Fj^+g—=0, 
OV ax OV OV 

which is the equation defining the Jacobi multipliers, because div X = dvF. 

Then, the inverse problem reduces to find the function g which is a Jacobi 
multiplier and L is obtained by integrating the function g two times with respect 
to velocities. The function L so obtained is unique up to addition of a gauge term. 

Coming back to the second-order Abel equation case, one can easily check that 
T>i{x, v) = V + x^ is a. Darboux polynomial for F with cofactor — since 

V ^ + Fa ^ ) {v + x^) = -x'^{v + x^) . 
Ox OV J 

The divergence of the vector field F is —Ax^, and then we see that there is a 
Jacobi multiplier of the form R = 2?^'*. Consequently, the Abel equation admits 
a Lagrangian description by means of a function L such that 

d^L 3 4 

— — = l^v + X ) , 

from where we obtain the Lagrangian L = La given by (|17p . 

But T>2{x, v) = 'iv + x^ is also a Darboux polynomial for F with cofactor — 3x^, 

V + Fa (3v + x^) = 3x^v — S(4:X^v + x^) = —Sx^(3v + x^) , 
Ox OV J 

and then we can find another Jacobi multiplier of the form I?2^ with ^2 = —4/3. 
The Abel equation admits a Lagrangian description by means of a function L such 
that 

d L 3-,-4/3 

9^ = (3- + -) , 
from where we obtain the Lagrangian L = La given by (|18p . 
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